Phase field model of premelting of grain boundaries 
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We present a phase field model of solidification which includes the effects of the crystalline orien- 
tation in the solid phase. This model describes grain boundaries as well as solid-liquid boundaries 
within a unified framework. With an appropriate choice of coupling of the phase field variable to the 
gradient of the crystalline orientation variable in the free energy, we find that high angle boundaries 
undergo a premelting transition. As the melting temperature is approached from below, low angle 
grain boundaries remain narrow. The width of the liquid layer at high angle grain boundaries di- 
verges logarithmically. In addition, for some choices of model coupling, there may be a discontinuous 
jump in the width of the fluid layer as function of temperature. 
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I. INTRODUCTION 

Understanding the structure of a boundary between 
two crystalline grains is a prerequisite to building theo- 
ries of microstructure development and evolution. Even 
in pure materials, there are many possibilities for this 
structure. Whereas a low angle grain boundary may be 
thought of as an array of dislocations , this description 
is not always accurate or useful. Several molecular dy- 
namics studies indicated that grain boundaries develop 
a layer of disordered, amorphous material as the tem- 
perature is raised || ^ H. Within a few degrees of 
the melting temperature Tm, grain boundaries in pure 
aluminum liquefy j^. The directly measured width of 
the liquid layer was consistent with a divergence at T^. 
Careful calorimetry measurements as in j?) may be able 
to pin down the nature of the divergence. Measurements 
of the dihedral angle of the bicrystal in contact with the 
melt in bismuth |^, ^ uncovered a discontinuous tran- 
sition as a function of misorientation of tilt boundaries 
near the melting point. 

Indirect probes of the structure of grain boundaries 
also find evidence of a structural transition. Measure- 
ments of grain boundary mobility pT| |, shear resistance 
and diffusion coefficients all found discontinuous 
jumps in these quantities at some characteristic temper- 
ature below the melting point. 

Melting at a grain boundary as well as surface melting 
[Q is a particular case of a broad class of wetting phe- 
nomena [y_5| . Generally, interfacial melting may be com- 
plete (the width of the liquid layer diverges as T ^ T^) 
or incomplete (it stays finite) depending on the functional 
dependence of the free energy on the thickness of the liq- 
uid layer. There may also be a discontinuous jump in the 
width of the liquid layer as a function of temperature. 

An explicit calculation of the free energy of a liquid 
layer is difficult and has been performed only in a few 
cases. For example, Kikuchi and Cahn calculated the free 
energy of a disordered grain boundary within a lattice gas 
model . They concluded that when the interaction of 
the solid-liquid interfaces is short range, the transition. 



if it exists, is continuous and the width of the liquid film 
diverges as a logarithm of the undercooling. Elbaum and 
Schick calculated the energy of a water film on ice due to 
the van der Waals dispersion forces and concluded 
that premelting there is incomplete. 

In this paper we study the modification of the phase 
field model of grain boundaries ||l8|. Our model treats 
solidification as well as grain boundaries within a unified 
consistent framework. Therefore, we ought to be able to 
describe melting of grain boundaries. The present au- 
thors elucidated the general properties of this model in 
Ref. There we were able to compute, in the limit 

of an infinitely sharp interface, the energy, width and 
mobility of model grain boundaries. 

We hope to use the predictions of our phase field model 
near the melting point to pin down the correct/realistic 
form of the model coupling functions which offer the es- 
sential flexibility in our model. We apply the results 
of @ and indeed find that grain boundaries undergo 
complete or incomplete continuous premelting transition 
depending on the misorientation of the two grains. In 
addition, for a particular choice of the model coupling 
functions, there is a discontinuous jump in the boundary 
width as a function of temperature. 



In the following Section II we briefly describe the phase 
field model of grain boundaries (for a better discussion 
see Ref. Q). In Sec. Ill we apply the results of Q 
to the modified model. We consider the structure com- 



posed of the grain boundary plus the layer of liquid and 
calculate its energy and width. We find that high angle 
grain boundaries undergo a continuous premelting tran- 
sition. We also find that for a particular choice of the 
model coupling functions, there is a discontinuous jump 
in the width of the liquid layer for some boundaries. In 
the next Sec. IV, we show that the existence of a liquid 
layer greatly affects the grain boundary migration. We 
summarize our findings in Sec. 
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II. MODEL 



The phase field is introduced in standard treatments 
of the kinetics of the Uquid-sohd transition ||2^ . It mea- 
sures the degree of structural disorder. Conventionally 
(j) — 1 corresponds to a perfect solid while cj) = corre- 
sponds to a liquid. In |Q Kobayashi Warren and Carter 
introduced rj which measures the degree of orientational 
disorder. We postulate that these two order parameters 
are essentially the same. Therefore the bulk free energy 
density /(0) in the free energy T 



T[<l), 9]^- [ dA 
e Jn 



/(0) + yIV0P 



-gWs\ve\ 

^2 



(1) 



should be thought of as a double well with minima at 
(f) — 0, 1. Here 9 is the local crystal orientation and a, 
e and s are positive model parameters. The model cou- 
pling functions g{4>) and h{(j)) must vanish in the fluid 
(j) = since the orientation 9 has no meaning there and 
therefore its gradients should not be penalized. 

The present authors analyzed the gradient flow of this 
free energy 



dt 



a^V^(f>- U-g^s\V9\ 



he^W9 + gs 



V9_ 



(2a) 
(2b) 



in Ref. ||T^. We used the subscripts to denote differ- 
entiation except for the time constants and Tg. The 
mobility functions P and Q must be positive deflnite, but 
are otherwise unrestricted. 





FIG. 1: Order parameters in a flat stationary grain boundary 
as a function of the coordinate normal to the boundary. The 
orientation 6 varies smoothly between and Ad in the grain 
boundary ( € [—£,£] and is constant in the interior of the 
grains. The phase field (p approaches 1 exponentially away 
from the boundary. 



In Ref. |19| we computed the structure of the boundary, 
its energy and mobility the limit of a sharp interface. Let 
us briefly summarize these results here. The structure of 
the boundary is shown in Fig. ^ The values of the phase 
field (j) in the center of the grain boundary 4>min and at 
the edge of the grain boundary 0niax are found by solving 
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FIG. 2: Amount of liquid W in units of ti a function of tem- 
perature L for a low-angle grain boundary, h = (jp' , s = 1 and 
M/a = O.L 



the following two equations |19| 



g{(f>miLx) = g(0min) + 

A9 
25l 
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d(f) {g{(f)„ 



9) 



0mi„ h^2f - P (g((?!)max) - g)'^/h 



(3) 
, (4) 



where a — a/e and s — s/e. Once ^min and ^max are de- 
termined, we can compute the energy of the grain bound- 
ary 7gb and the amount of "fiuid" or disordered material 
at the interface W via 
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III. FLAT STATIONARY BOUNDARY 

Equipped with the way of calculating measurable prop- 
erties of the grain boundary in our model we make choices 
for the coupling functions and explore the results of the 
calculation. As we alluded to before the bulk free energy 
density /(</>) is a double well which we choose to be pO] 



/(0) = 



02(1 _^)2 



L{1 



5(1 + 30 + 602), (7) 



where L = — /(O) is linear in the temperature devia- 
tion from the melting point T — Tm. The choices of 
g{(f>) and h{(j)) are more flexible. These functions provide 
the essential freedom in our model. The only constraint 
on them is that they must vanish sufficiently quickly at 
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= 0. For example, they can be chosen [g8| to match 
the Read-Shockley [2^ misorientation dependence of the 
grain boundary energy. In this paper we wiU set 

5(0) = 0', (8) 

and explore two choices for h 

hi{(j)) = 02 and /i2(0) = 1. (9) 

The second choice is problematic since it ascribes an en- 
ergy cost to changes of orientation in the liquid where it 
has no meaning. We nevertheless explore its implications 
to gain insight and intuition into the role of the model 
function h. 
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FIG. 3: Energy 7gb in units of S of the low-angle grain bound- 
ary (same values of the parameters as in Fig. 0) and the saddle 
point solution (upper branch). 



A. h^(j>^ 

A typical dependence of the amount of liquid W on 
temperature L for a low-angle boundary is shown in 
Fig. 1^. The meaning of the "low-angle" will become 
apparent in a moment. Below the melting temperature 
L < 0, the narrow, dry boundary is the only solution. 
Above the melting temperature, the narrow, dry bound- 
ary is metastable with respect to the uniform fluid. The 
second solution (upper branch) is the saddle point on 
the path from the narrow boundary to the globally sta- 
ble uniform fluid. It gives the size of the critical nucleus. 
The energy barrier for nucleation ^7 shown in Fig. ^ van- 
ishes at the spinodal temperature Lg above which the 
grain boundary cannot exist. Note that since the width 
of the saddle point solution (critical nucleus size) diverges 
(logarithmically |2^) at the melting temperature, the en- 
ergy of the saddle point solution must approach that of 
the two liquid solid interfaces 

271^ = 25^ d0y27=|, (10) 

which it indeed does. 

The spinodal temperature's dependence on s and the 
misorientation AO is shown in Fig. 0. As expected the 
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FIG. 4: The spinodal temperature La, above which the grain 
boundary cannot exist, as a function of the misorientation. 
Note that the spinodal temperature vanishes at a critical mis- 
orientation. 



spinodal temperature vanishes at an s-dependent critical 
misorientation AOc- Therefore, high-angle grain bound- 
aries {A9 > A9c) cannot exist above the melting tem- 
perature. For these boundaries, the amount of fluid W 
diverges as the melting point is approached from below. 
Another way to define this critical misorientation is to 
notice that the energy 7gb of all high-angle grain bound- 
aries approach that of two liquid-solid interfaces ( |lO| ) at 
the melting point, whereas low-angle grain boundaries 
are less energetically costly. 
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FIG. 5: Critical misorientation in units of a as function of s. 
h = (fp' . Solid line is the analytic result in the limit of large s. 

The dependence of the critical misorientation on s, 
shown in Fig. ^, can be calculated approximately in the 
limit of large s. We need to find the largest misorien- 
tation for which a solution to (^) exists at the melting 
point L = Q. Using the approximate expression for the 
right hand side of Eq. (^, calculated in ||l9|, we obtain 
(0 subscript denotes differentiation) 

AO 2 77(0 

a S 5'</,(0min) S 

And since is restricted to lie within the [0, 1] interval. 
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a solution only exists for 



M < Mc = -. 

s 



(12) 



Indeed, for large s the numerically found critical misori- 
entation shown in Fig. ^ approaches this prediction. It 
is worth noting that no solutions to our model exist for 
s < 0.5. In this case, some values of 0min < 1 result in 
unphysical (/)inax according to Eq. 





1 1 

ASi = 0.655 




A02 = 0.5a 




/ ; 


A6»3 = 0.35a 














1 1 




1 1 



-0.02 



-0.01 



0.01 



0.02 



FIG. 6: Amount of liquid W for h = \va units of a as function 
of temperature L for s = 1 and three different misorientations 
A6l„ < ASi, AS, < Ae2 < A6l„ and A6I3 < Ae,. 



B. 



h = 1 



The behavior of W is more complex in this case. As 
shown in Fig. ^ there are three distinct cases. The 
amount of liquid at a high-angle boundary diverges 
smoothly at the melting point. When is below an 
upper critical misorientation A^^, there is a range of tem- 
peratures below melting in which there are two solutions, 
one stable and one metastable. For example, by looking 
at the energies of these solutions, we discover, that the 
narrow, dry solution is stable up to the transition temper- 
ature Lt < (see Fig. Above this transition temper- 
ature, the wide, wet solution is stable, while the narrow, 
dry solution is metastable. The transition from narrow, 
dry boundary to wide, wet boundary is discontinuous. 

Below a lower critical misorientation A9i , there is still 
a region of temperatures in which there exists a wide, 
wet solution, but the narrow, dry solution is stable all 
the way up to the melting temperature. 

This definition of the upper and lower critical misorien- 
tations lends itself well to predicting the result of raising 
the temperature slowly. According to our definition, a 
low-angle grain boundary (A6 < A0i) will remain nar- 
row up the melting point where it will melt via a first 
order activated transition. Boundaries with intermedi- 
ate misorientations (A^; < A9 < A0„) will experience a 
discontinuous jump in width at some characteristic tem- 
perature Lt < 0. When the temperature is increased 




FIG. 7: Upper and lower critical misorientations as a function 
of s for the h = 1 case. 



further, the width of the liquid layer will diverge loga- 
rithmically at the melting point. High- angle boundaries 
{A9 > A9u) experience a smooth divergence of the width 
of the liquid layer as the melting point is approached. 
The dependence of the lower and upper critical misori- 
entations on s is shown in Fig. 0. 



IV. GRAIN BOUNDARY MOBILITY 

In this paper we would like to illustrate the flexibility 
of our model by demonstrating that the choice of model 
functions is sufficient to reproduce a wide range of be- 
haviors observed experimentally and numerically. For 
example, molecular dynamics simulations as well as 
experiments reveal that the mobility of tilt grain 
boundaries vanishes quickly at the zero misorientation. 
The same studies find (in two dimensions) a sharp peak 
in the mobility at special (high coincident site density) 
misorientations at which the energy of the boundary has 
a cusp. 

As we showed in Ref. within our model, curved 
boundaries move with a normal velocity proportional to 
their curvature, energy 7gb and mobility M. which we 
computed in the sharp interface limit 
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/lV2/-S2(5(0max)-5)V/^ 

gV2/-S2(g(</,,„ax)~g)2A 



a 



(13) 



where f ^ = / and fg — Tg/e^. 

We studied the scaling of the mobility in the limit of 
the vanishing misorientation A9 and found that this limit 
is governed by the behavior of the model functions near 
(j) = 1. In this paper g and h are regular at </> = 1. There- 
fore, the behavior of P and Q near 0=1 determines the 
scaling of the grain boundary mobility M for small mis- 
orientations. If the most strongly divergent one of the 
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two functions P and Q behaves like (1 — (t>)~^ , then the 
mobiUty scales like (A0)^^^. As we remarked in the 
divergence of P and Q at = 1 is physically plausible 
since we expect the dynamics of (j) and 9 to slow down in 
the grain interior. 
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FIG. 9: Mobilities of the wet and dry boundary solutions. 
Here s — & = 1, AS — 0.55 and h = 0^. 
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FIG. 8: Grain boundary mobility M as, a. function of misori- 
entation both in units of a for h — (p^ and three choices of P. 
The slopes of the two straight line segments are +1 and —1. 
Temperature is L = —0.001, s = a = 1. Note that above the 
critical angle, different P makes does not influence the grain 
boundary mobility. 

In Fig. ^ we show the grain boundary mobility as a 
function of misorientation ior h — (p^ , Q — — log(l — (p) 
and three different choices of P. Indeed, the behavior of 
P near 0=1 dictates whether the mobility Ai vanishes 
or diverges at zero misorientation. The undercooling L = 
—0.001 is small in this figure. Larger undercoolings do 
not alter the qualitative features of the situation. Our 
model can therefore be adjusted to describe systems in 
which the mobility dips at small misorientations as well 
as systems which exhibit peaks in the mobility at special 
misorientations . 

Since the grain boundary mobility is strongly depen- 
dent on the amount of disorder W, its behavior as a func- 
tion of temperature can be inferred from its behavior as 
a function of misorientation. Specifically, when /i > 2, 
the mobility is an increasing function of the amount of 
disorder. And vice versa is a decreasing function of 
W when fj, < 2. The amount of disordered material at 
the grain boundary grows with increasing temperature 
as well as with increasing misorientation. We therefore 
expect the mobility to be an increasing function of tem- 
perature when fi > 2. 

To illustrate this point we set h = 1 and plot in Fig. ^ 
the mobility of a grain boundary which undergoes a dis- 
continuous premelting transition, corresponding to A^2 
in Fig. I {Adi < M2 < A6'„). We also choose /i = 3 
so that narrow boundaries move slower. The mobility of 
the narrow, dry boundary, which is stable below Lj, is 
significantly smaller than that of the wide, wet bound- 
ary. Thus when the boundary undergoes the transition 
it will suddenly accelerate to almost twice the speed. 



V. CONCLUSIONS 

In this article we explored the properties of a phase 
field model of crystal grains. The model is constructed 
by introducing an order parameter (j) which measures the 
degree of disorder (or "fiuidity") and the average local 
orientation 9. When the phenomenological free energy 
contains a term which is linear in |V0|, its relaxational 
kinetics is singular and has to be interpreted with the 
help of the extended gradient theory. The relaxing so- 
lutions thus obtained are interpreted as a collection of 
grains in which 9 is spatially uniform, separated by nar- 
row grain boundaries. 

Our main goal was to demonstrate the ability of this 
phase field model to reproduce experimentally observed 
behaviors of grain boundaries. Since our model treats 
solid grain boundaries and liquid-solid interfaces within 
a unified framework, it ought to correctly describe the 
premelting transition of grain boundaries. We indeed 
find that depending on the choice of the model functions 
the grain boundary undergoes either a continuous or a 
discontinuous premelting transition. 

If we choose the coefficient of (V0)^ in the free en- 
ergy to be h{(j)) = (jP, the transition is continuous. This 
means that there exists a critical misorientation A9c such 
that grain boundaries with lower misorientations remain 
narrow and dry all the way up to the melting temper- 
ature. When the misorientation is larger than critical, 
grain boundaries develop a layer of fluid whose width di- 
verges at the melting temperature as the logarithm of 
the undercooling. The logarithmic divergence is a con- 
sequence of the short range interaction of liquid-solid in- 
terfaces within this theory. 

When h — 1, there is a range of misorientations for 
which the width of the liquid film undergoes a jump at 
some characteristic temperature below the melting tem- 
perature. If the temperature is raised further, the width 
of the liquid film diverges logarithmically as before. Low 
angle grain boundaries remain dry and narrow as before. 
The jump in the width of the fiuid layer is due to the de- 
tails of the short-range interaction of the two liquid solid 
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interfaces and should be observable in some real systems. 

Finally, we have shown that by adjusting the nature 
of the divergence of the model function P we are able to 
obtain a wide range of behaviors of the grain boundary 
mobility. In particular, depending on the exponent of 
the divergence of P near 0=1, the mobility of the grain 
boundary can either vanish at the zero misorientation or 
it can diverge. Thus the model can emulate the behavior 
of the mobility observed experimentally near the special 



misorientations as well as near the zero misorientation. 

In principle, by introducing appropriate additional 
terms in the free energy, one can construct a more re- 
alistic model in which the grain boundary energy has 
many minima at special misorientations in addition to 
the cusp at the zero misorientation. It would then be 
possible to choose P and Q to ensure appropriate peaks 
or zeroes of the grain boundary mobility at the special 
misorientations as well as the zero misorientation. 
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